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ABSTRACT: In this paper we introduced an approach to increase integration rate of element of a 


conventional energy efficient double tail dynamic comparator. We formulate recommendations for 
optimization of annealing to decrease dimensions of field-effect heterotransistors and to increase their 
density. At the same time one can decrease transition area between doped and undoped areas of 
elements of the considered comparator. We introduce an analytical approach to model diffusion and 
ion types of doping with account concurrent changing of parameters in space and time. At the same 
time the approach gives us possibility to take into account nonlinearity of considered processes. 
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INTRODUCTION 

Actual and intensively solving problem in 
solid state electronics is to increase the 
integration rate of elements of the integrated 
circuits (p-n-junctions, their systems and so on) 
[1-8]. Increasing the integration rate leads 
necessity to decrease their dimensions. Several 
approaches were adopted to decrease the 
dimensions. They are generally using laser and 
microwave types of annealing of infused 
dopants. These types of annealing are also 
widely using for annealing of radiation defects, 
generated during ion implantation [9-17]. Using 
the approaches gives a possibility to increase 
integration rate of elements of integrated 
circuits through inhomogeneity of technological 
parameters due to generating inhomogenous 
distribution of temperature. In this situation 
one can obtain decreasing dimensions of 
elements of integrated circuits [18] with 
account Arrhenius law [1,3]. Another approach 
is to manufacture elements of integrated 
circuits with smaller dimensions is doping of 
heterostructure by diffusion or ion implantation 
[1-3]. However in this case optimization of 
dopant and/or radiation defects is required 
[18]. 

In this paper, we consider a heterostructure. 
The heterostructure consist of a substrate and 


several epitaxial layers. Some sections have 
been manufactured in the epitaxial layers. 
Further we consider the doping of these 
sections by diffusion or ion implantation. The 
doping gives a possibility to manufacture field- 
effect transistors framework a conventional 
energy efficient double tail dynamic comparator 
so as it is shown on Figs. 1. The manufacturing 
gives a possibility to increase density of 
elements of a conventional energy efficient 
double tail dynamic comparator [4]. After 
considered doping dopant and/or radiation 
defects should be annealed. Framework the 
paper we analyzed dynamics of redistribution 
of dopant and/or radiation defects during their 
annealing. We introduce an approach to 
decrease dimensions of the element. However it 


is necessary to complicate technological 
process. 
METHOD OF SOLUTION 


In this section we determine spatio-temporal 
distributions of concentrations of infused and 
implanted dopants. To determine these 
distributions we calculate appropriate solutions 
of the second Fick's law [1,3,18] 
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Figure 1 The considered the comparator [4] 
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The function C(x,y,z,t) describes the spatio- 
temporal distribution of concentration of 
dopant; T is the temperature of annealing; Dc is 
the dopant diffusion coefficient. Value of dopant 
diffusion coefficient could be changed with 
changing materials of heterostructure, with 


changing temperature of materials (including 
annealing), with changing concentrations of 
dopant and radiation defects. We approximate 
dependences of dopant diffusion coefficient on 
parameters by the following relation with 
account results in Refs. [20-22] 
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Here the function D; (x,y,z,T) describes the 
spatial (in heterostructure) and temperature 
(due to Arrhenius law) dependences of diffusion 
coefficient of dopant. The function P (x,y,z,T) 
describes the limit of solubility of dopant. 
Parameter vy €[1,3] describes average quantity 
of charged defects interacted with atom of 
dopant [20]. The function V (x,,z,t) describes 
the spatio-temporal distribution of 
concentration of radiation vacancies. Parameter 


"(ay20) 
P’(x,y,z,T) 


E-ISSN: 2349 5359; P-ISSN: 2454-9967 


V"(x, y,z,t) 
V S> vy 


V’ describes the equilibrium distribution of 
concentration of vacancies. The considered 
concentration dependence of dopant diffusion 
coefficient has been described in details in [20]. 
It should be noted that using diffusion type of 
doping did not generation radiation defects. In 
this situation, G7= ¢2= 0. We determine spatio- 
temporal distributions of concentrations of 
radiation defects by solving the following 
system of equations [21,22] 


(3) 


Ot Ox 


Ox 


OI (x, y,zt)_ 0 o,G.y.2.7) eee), 2 b,(sy.27) Hee), 
y 


Oy 


~ a cs yor ened) —k,,(x,y,z,T) U(x, y,z,t) V(x, y,z,t)- 


z z 
—k, (x, y,z,T I(x, y,z,t) (4) 
OV (x,y,z,t)_ 2 OV(x,y,z,t)|, 2 OV (x, y,z,t) 
= Dx, y,z,T D,\x, y, z,T )—-————- 
os 2 v(% 95.7) . (x, y,z,7) o- 


+2 vloyer) VELA |, ies 2 VGoe: 
xz j 


+k, (~%y,zTV"(~y,z0). 


Boundary and initial conditions for these equations are 
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Here p =LV. The function I (x,y,z,t) describes the 
spatio-temporal distribution of concentration of 
radiation interstitials; D,(x,,z,T) are the 
diffusion coefficients of point radiation defects; 
terms V2(x,y,z,t) and /[2(x,y,z,t) correspond to 
generation divacancies and _ diinterstitials; 
kiyv(X%y,z,T) is the parameter of recombination of 
point radiation defects; ki(xy,z,T) and 


kvv(xy,z,T) are the parameters of generation of 
simplest complexes of point radiation defects. 


Further we determine distributions in space 
and time of concentrations of divacancies 
O/(x,y,z,t) and diinterstitials @(x,y,z,t) by 
solving the following system of equations 
[21,22] 
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Boundary and initial conditions for these equations are 
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D (x,y,2,0)=far (x,y,z), Dy (x,y,2,0)=fov (xy,2Z). (7) 


Here Do,(x,y,z,T) are the diffusion coefficients of 
the above complexes of radiation defects; 
ki(xy,z,T) and ky (x,y,z,T) are the parameters of 
decay of these complexes. 


We calculate distributions of concentrations 
of point radiation defects in space and time by 
recently elaborated approach [18]. The 
approach based on_ transformation of 
approximations of diffusion coefficients in the 
following form: D,(x,,z,T)=Dop[1+ & gp{xy,Z,T)], 
where Do, are the average values of diffusion 


I(xy,zt)=Iayztyl, V(xy,zt)= 


coefficients, 0<é,<1, |g,(x, y,z,T)|<1, e =LV. We 
also. used analogous transformation of 
approximations of parameters of recombination 
of point defects and parameters of generation of 
their complexes: kiv(Xy,Z,T)=Korv[1+é,v 
Giv(Xy.Z,T)], ki(Xy,2,T)=ko [1+e1 gui(%y,2,T)] 
and kyv(x,y,z,T)=Kovv [1+ evv gv.v(%y,z,T)], where 
kopi2 are the their average values, O<e,y <1, 
O<e, <1, Osevv<l1, | givlxy,z,T)|<1, | 
gi(Xy,Z,T)|<1, \gvv(Xy,zZ,T)|<1. Let us introduce 
the following dimensionless variables: 
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transformation of Eqs.(4) and conditions (5) to the following form 
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We determine solutions of Eqs.(8) with conditions (9) framework recently introduced approach [18], 
i.e. as the power series 


P(1.0.9,9)= Le, Lo! LOD, (11.6.9). (10) 
i= j= = 
Substitution of the series (10) into Eqs.(8) and Now we can calculate the distributions of 
conditions (9) gives us possibility to obtain concentrations of simplest complexes of point 
equations for initial-order approximations of radiation defects in space and time. To 


determine the distributions we transform 

approximations of diffusion coefficients in the 

(wn, d, 9) and corrections for them following form: Do(Xy,Z,T)=Doe@[1+ 
E09 0(XY,Z,1)], where Dog are the average 

Ly 1); p, 3) and V,.(757,99), i21,j21,k Scene a soetficlents In this ann 

>1. The equations are presented in the the Eqs.(6) could be written as 

Appendix. Solutions of the equations could be 

obtained by standard Fourier approach [24,25]. 

The solutions are presented in the Appendix. 
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Further we determine solutions of above equations as the following power series 


® (x, y,z,t)= ve 


Now we have used the series (11) into Eqs.(6) 
and appropriate boundary and _é initial 
conditions. The using gives the possibility to 
obtain equations for initial-order 
approximations of concentrations of complexes 
of defects @,0(x,y,z,t), corrections for them 
@,(x,y,z,t) (for them i >1) and boundary and 
initial conditions for them. We _ remove 
equations and conditions to the Appendix. 
Solutions of the equations have been calculated 
by standard approaches [24,25] and presented 
in the Appendix. 


® ,,(x,y,z,t). (11) 


Now we calculate the distribution of 
concentration of dopant in space and time by 
using the approach, which was used for analysis 
of radiation defects. To use the approach we 
consider following transformation of 
approximation of dopant diffusion coefficient: 
D,(%y,Z,T)=Dor[1+ é1gi(X,y,Z,T)], where Doz is the 
average value of dopant diffusion coefficient, 
O<ex< 1, |gi(xy,z,T)|<1. Farther we consider 
solution of Eq.(1) as the following series: 


C(x, y, z,t)= were Cs V; z,t) 


Using the relation into Eq.(1) and conditions (2) 
leads to obtaining equations for the functions 
Ci(xy,z,t) (i 21, 7 21), boundary and initial 
conditions for them. The equations are 
presented in the Appendix. Solutions of the 
equations have been calculated by standard 
approaches (see, for example, [24,25]). The 
solutions are presented in the Appendix. 


We analyzed distributions of concentrations 
of dopant and radiation defects in space and 
time analytically by using the second-order 
approximations on all parameters, which have 
been used in appropriate series. Usually the 
second-order approximations are enough good 
approximations to make qualitative analysis 
and to obtain quantitative results. All analytical 
results have been checked by numerical 
simulation. 


DISCUSSION 


In this section we analyzed spatio-temporal 
distributions of concentrations of dopants. Fig. 


concentrations of dopants in neighborhood of 
interfaces of heterostructures. We calculate 
these distributions of concentrations of dopants 
under the following condition: value of dopant 
diffusion coefficient in doped area is larger, than 
value of dopant diffusion coefficient in nearest 
areas. In this situation one can find increasing of 
compactness of field-effect transistors with 
increasing of homogeneity of distribution of 
concentration of dopant at one time. Changing 
relation between values of dopant diffusion 
coefficients leads to opposite result (see Figs. 3). 


It should be noted, that framework the 
considered approach one_ shall optimize 
annealing of dopant and/or radiation defects. 
To do the optimization we used recently 
introduced criterion [26-34]. The optimization 
based on approximation real distribution by 
step-wise function yw (x,y,z) (see Figs. 4). Farther 
the required values of optimal annealing time 
have been calculated by minimization the 
following mean-squared error 


2 show _ typical sa distributions of 
X,Y,Z; x,y, z)|dzdydx. 12 
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Fig. 2a. Dependences of concentration of dopant, infused in heterostructure from Figs. 1, on 
coordinate in direction, which is perpendicular to interface between epitaxial layer substrate. 
Difference between the values of dopant diffusion coefficient in layers of heterostructure increases 
with increasing of number of curves. Value of dopant diffusion coefficient in the epitaxial layer is 
larger, than value of dopant diffusion coefficient in the substrate 


2.0 


T T — 
0 L/4 L/2 3L/4 L 
x 


Fig. 2b. Dependences of concentration of dopant, implanted in heterostructure from Figs. 1, on 
coordinate in direction, which is perpendicular to interface between epitaxial layer substrate. 
Difference between values of dopant diffusion coefficient in layers of heterostructure increases with 
increasing of number of curves. Value of dopant diffusion coefficient in the epitaxial layer is larger, 
than value of dopant diffusion coefficient in the substrate. Curve 1 corresponds to homogenous sample 
and annealing time © = 0.0048 (L,2+L,2+L,2)/Do. Curve 2 corresponds to homogenous sample and 
annealing time © = 0.0057 (L,2+Ly2+L,2)/Do. Curves 3 and 4 correspond to heterostructure from Figs. 
1; annealing times © = 0.0048 (L,2+Ly2+L,2)/Do and © = 0.0057 (L,* +Ly?+L,2)/Do, respectively 
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Fig.3a. Distributions of concentration of dopant, infused in average section of epitaxial layer of 
heterostructure from Figs. 1 in direction parallel to interface between epitaxial layer and substrate of 
heterostructure. Difference between values of dopant diffusion coefficients increases with increasing 
of number of curves. Value of dopant diffusion coefficient in this section is smaller, than value of 
dopant diffusion coefficient in nearest sections 


f(x) 


1.00000 


0.10000 


Epitaxial Substrate 


0.01000 layer 2 


0.00100 


C(x,0) 


0.00010 


0.00001 


' Epitaxial 
| layer 1 


0.00000 T T T 
0 L/4 LI2* 3L/4 L 
x 


Fig.3b. Calculated distributions of implanted dopant in epitaxial layers of heterostructure. Solid lines 
are spatial distributions of implanted dopant in system of two epitaxial layers. Dushed lines are spatial 
distributions of implanted dopant in one epitaxial layer. Annealing time increases with increasing of 
number of curves 
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Fig.4a. Distributions of concentration of infused dopant in depth of heterostructure from Fig. 1 for 
different values of annealing time (curves 2-4) and idealized step-wise approximation (curve 1). 
Increasing of number of curve corresponds to increasing of annealing time 
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Fig.4b. Distributions of concentration of implanted dopant in depth of heterostructure from Fig. 1 for 
different values of annealing time (curves 2-4) and idealized step-wise approximation (curve 1). 
Increasing of number of curve corresponds to increasing of annealing time 
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Fig. 5a. Dimensionless optimal annealing time of infused dopant as a function of several parameters. 
Curve 1 describes the dependence of the annealing time on the relation a/L and € = y= 0 for equal to 
each other values of dopant diffusion coefficient in all parts of heterostructure. Curve 2 describes the 
dependence of the annealing time on value of parameter ¢ for a/L=1/2 and €= y= 0. Curve 3 describes 
the dependence of the annealing time on value of parameter & for a/L=1/2 and ¢ = y= 0. Curve 4 
describes the dependence of the annealing time on value of parameter yfor a/L=1/2 and e= €=0 


We show optimal values of annealing time as 
functions of parameters on Figs. 5. It is known 
that standard step of manufactured ion-doped 
structures is annealing of radiation defects. In 
the ideal case after finishing the annealing 
dopant achieves interface between layers of 
heterostructure. If the dopant has no enough 
time to achieve the interface, it is practicably to 

0.12 


0.08 


@©D,L’ 


0.00 


anneal the dopant additionally. The Fig. 5b 
shows the described dependences of optimal 
values of additional annealing time for the same 
parameters as for Fig. 5a. Necessity to anneal 
radiation defects leads to smaller values of 
optimal annealing of implanted dopant in 
comparison with optimal annealing time of 
infused dopant. 
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Fig.5b. Dimensionless optimal annealing time of implanted dopant as a function of several parameters. 
Curve 1 describes the dependence of the annealing time on the relation a/L and € = y= 0 for equal to 
each other values of dopant diffusion coefficient in all parts of heterostructure. Curve 2 describes the 
dependence of the annealing time on value of parameter ¢ for a/L=1/2 and €= y= 0. Curve 3 describes 
the dependence of the annealing time on value of parameter & for a/L=1/2 and ¢ = y= 0. Curve 4 
describes the dependence of the annealing time on value of parameter yfor a/L=1/2 and e= €=0 


CONCLUSION heterotransistors and to increase their density. 
At the same time one can decrease transition 
area between doped and undoped areas of 
elements of the considered comporator. We 
introduce an analytical approach to model 
diffusion and ion types of doping with account 


concurrent changing of parameters in space and 


In this paper we introduced an approach to 
increase integration rate of element of a 
conventional energy efficient double _ tail 
dynamic comparator. We formulate 
recommendations for optimization of annealing 
to decrease dimensions of field-effect 
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time. At the same time the approach gives us 
possibility to take into account nonlinearity of 
considered processes. 
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Equations for functions ®,;(x,y,z,t), 120 to describe concentrations of simplest complexes of radiation 
defects. 
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t Ox Oy OZ 
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PA g, (x,y,z, OL X-2) = O g, (x. y,2,7 Lule r24) , 
Oy Oy Oz Az 
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Coo(x,y,2,0)=fc (xy,Z), Ci(xy,Z,0)=0, i 21, j =1. 


Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.9 No.1 2517-2538 (2022) 2536 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


Functions C;(x,y,z,t) (i 20, j 0) could be approximated by the following series during solutions of the 
above equations 


Fie — 
Cuo(X ¥ Zt) = LLL LLL: DY Fc (x)e, (ye, (zenclt): 
Hereeé, A(t) = exp| — 1°n’D, ae Al 
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dwdvdudt; 
mt V,W, T) Ow 
Cie me)= ppg Enhaealel Pee eee 


dwdvdud 


)2Cu(uv.w.r) 


x u,Vv,w,l 
g,( Ou 


Li (x)e,(y)e,(Z)eyclt) x 


x Je,c(- r)fe,(u) fs O)fer(we,(us.w.7)2S (aa — * 


* mn e (tie (- ele leo)fs,0v)g, (2,07) 9) a wavdude * 
0 0 W 


« Faceladey(y)e,(@)- gE Ft einai ele) fssle)fes(r) 


L, C’ (u,v, w,T) 6C,,(u,v,w, t) 
xn Jo) Gav wT) Ou ne “ ‘i 


Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.9 No.1 2517-2538 (2022) 2538 E-ISSN: 2349 5359; P-ISSN: 2454-9967 
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